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We study the phenomenon of migration of the small molecular weight component of a binary 
polymer mixture to the free surface using mean field and self-consistent field theories. By proposing a 
free energy functional that incorporates polymer-matrix elasticity explicitly, we compute the migrant 
volume fraction and show that it decreases significantly as the sample rigidity is increased. Estimated 
values of the bulk modulus suggest that the effect should be observable experimentally for rubber¬ 
like materials. This provides a simple way of controlling surface migration in polymer mixtures 
and can play an important role in industrial formulations, where surface migration often leads to 
decreased product functionality. 
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Introduction: When a polymer mixture having mobile 
components of different molecular weights and with an 
interface that is free to atmosphere is left to equilibrate, 
the small molecular weight component migrates to the 
surface m- Several industrial formulations, e.g. choco¬ 
late [1], food packaging etc. suffer from this ubiqui¬ 
tous problem. While many experimental m 116] and 
theoretical studies [7] of this phenomenon exist, a good 
quantitative agreement between theoretical predictions 
and experiments is still lacking [T]. Further, experimen¬ 
tal strategies to control the amount of material migrating 
to the surface is in a nascent stage of development. 

In this letter we ask: How does the elasticity of 
the polymer matrix influence surface migration of small 
molecules in polymer mixtures? We propose a free en¬ 
ergy functional that incorporates elasticity of the poly¬ 
mer mixture explicitly, a feature that has been ignored 
in previous surface segregation studies. Using mean field 
theory (MFT) and self-consistent field theory (SOFT) 
we show that as the sample rigidity is increased (i) the 
migrant fraction decreases, and (ii) a wetting transition 
can be avoided (demonstrated by a geometric construc¬ 
tion lii). These results are of paramount importance in 
industrial product formulations where surface migration 
of small molecular weight component results in decreased 
functional performance of the product. 

Surface Migration: For a binary mixture, the compo¬ 
nent with the lower surface energy will migrate to the 
interface. A balance between loss of translational en¬ 
tropy and gain in surface energy dictates the equilibrium 
morphology of such systems. This is shown in Fig. [^with 
a high migrant (black) concentration close to the inter¬ 
face {z = 0) of a mixture of low and high (red) molecular 
weight polymers. The migrant concentration decreases 
monotonically to the bulk concentration (fao as z ^ oo. 
A crucial parameter that dictates the thermodynamics 
of the system is yiV, where x is the miscibility parame¬ 


ter, and N the molecular weight of the migrant. As 
increases a wetting transition, characterised by a macro¬ 
scopic wetting layer is observed (Fig. inset) [1]. 

Surface migration was first observed using X-ray pho¬ 
toemission spectroscopy nni and the resolution of the 
depth profile of the migrant concentration was improved 
significantly using neutron-reflectivity m- Further stud¬ 
ies concentrated on the theoretical aspects of the mi¬ 
gration by Schmidt and Binder m and subsequently 
a comparison between theory and experiments m- 
The wetting transition of polymer mixtures at the air- 
mixture interface was first demonstrated by Steiner et 
al. |14] . Experimental and theoretical developments of 
this phenomenon has recently been reviewed by a few 
authors HIS]. 

Flory-Huggins theory: The thermodynamics of mixing 
of two chemically different polymers is well described by 
Flory-Huggins (FH) theory [15]. The mixing free energy 
per unit volume is given by 

where y is the miscibility parameter, and Na, and Nb 
are the degree of polymerisation of A and B polymers 
respectively. The volume fractions of the A {(j)A = </>), 
and B {(j)B = ^ polymers in Eq. [^thus satisfy the 
incompressibility constraint (j)A + 4>b = 1 - The phase 
behaviour of such systems is well known m- Below 
a critical value of the miscibility parameter y < yc = 
1/(2A’^)-|-1/(2 A’b)-|-1/(-\/A’aA^b) the equilibrium phase 
is a homogeneous mixture of A and B polymers. For 
y > Xc {e.g. effected by changing temperature) phase 
segregation occurs with the system separating into A and 
B rich regions. Depending on the parameters, a first 
or second order transition might be observed. This is 
schematically shown in the inset of Fig. (solid line). 

Schmidt-Binder formalism: While FH free energy de¬ 
scribes the phase separation process in bulk it cannot be 
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directly applied to study segregation close to an interface 
that is exposed to atmosphere. Cahn’s [8] seminal work 
provides a cue in this case. This framework offers a way of 
calculating the concentration profile of a fluid near a wall, 
given a limiting fluid concentration, using the calculus of 
variations. The Flory-Huggins as well as Cahn’s theory 
have successfully been combined into a single mean field 
description to describe the surface segregation of binary 
polymer mixtures by Schmidt and Binder |12j (referred 
as SB henceforth). The SB free energy functional for a 
semi-infinite system (z > 0) with an attractive surface 
having area A at z = 0 is given by 
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where k{(j)) = is the coefficient associated 

with the energetic cost of creating a concentration gra¬ 
dient (obtained within the random phase approxima¬ 
tion [nma [16]), and A/i models the exchange chemi¬ 
cal potential. The SB functional also incorporates the 
surface free energy gain of the migrant Fs{4>i) expressed 
as a polynomial expansion of the migrant volume frac¬ 
tion at the surface, (c^i = 0(z = 0)) and is given by 
Fs{(j)i) = where is the surface chemical 

potential and the coefficient g characterises the change in 
bulk interactions due to the surface [T^[T7]. Within the 
gamut of square gradient theory the free energy func¬ 
tional in Eq.(|^ can be minimised SFsb[4']/^4' = 0, to 
yield an integral expression for z((j)), which can be in¬ 
verted to obtain the concentration profile of the migrant 
[12] . For small values of xfV an exponentially decay¬ 
ing profile shows reasonable agreement with experimental 
data [T]. 

Elastic Flory-Huggins theory: We now explore the role 
of polymer matrix elasticity in the small molecule migra¬ 
tion process. If one component (B polymer in our case) 
forms an elastic network as in cross-linked gels (retic¬ 
ulated permanent network) or transient network (as in 
polymer solutions) then its entropic contribution to the 
FH mixing free energy would be negligible in comparison 
to that of the migrant. Assuming Flory-Rehner form of 
free energy |18j describing the energy cost of a migrat¬ 
ing oligomer as it pushes its way through the matrix, the 
Flory Huggins elastic free energy Ffhe can be written as 
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FIG. 1. Schematic figure showing a mixture of low (black) and 
high (red) molecular weight polymers, with the low molecu¬ 
lar weight component migrating to the free interface 2 = 0. 
A semi-infinite geometry is assumed. The volume fraction of 
the migrant in bulk and at the surface is denoted by <j)oo, and 
4>i respectively. Inset shows migrant concentration profiles 
for different values of ■ For low values of yA a mono- 
tonically decreasing concentration profile is observed (dashed 
line). As yA increases a wetting transition, characterised by 
a macroscopically thick migrant layer (solid line with a break) 
is observed. 


4>b = VbIV respectively, we obtain an expression of 
the relative chain extension A in terms of the volume 
fractions, i.e. 4 >b/4>Bo = VbIVbo = [R/RaY = 
Since (j)Bo = 1 - '(’Aq = 1 - <^oo and (pB = 1-0, 
A = [(1 — 0)/(l — 0oo)]3. The volume fraction of mi¬ 
grant molecules in the sample is 0. Following the theory 
of elastic rubber networks m, the number of chains in 
the network can be estimated as n = Nb/V = (1 — 0oo)- 
The number of chains that form the polymer matrix The 
elastic free energy in Eq.(|^ is thus 
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where B is the elastic modulus expressed in units of fcsT. 
The free energy that describes the small molecule mi¬ 
gration through a matrix where elastic effects have been 
explicitly incorporated is therefore given by 
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where F^i = B^{X^ -I- | — 3), modelling uniaxial network 
deformation US], with A representing the relative chain 
extension, n the number of chains in the network and B 
the elastic modulus. The relative chain extension is de¬ 
fined as A = R/Ro, where Rq and R denote the length 
of the polymer before and after deformation respectively. 
Assuming the volume fraction of the matrix polymer be¬ 
fore and after deformation being pBo = VBo/Vy and 
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where Ffhe is the elastic Flory Huggins functional in 
Eq.(|^ and the gradient, exchange chemical potential and 
surface contributions to the free energy is the same as the 
SB free energy functional in Eq.([^. 

The role of elasticity in the phase separation of binary 
polymer mixtures where both species are cross-linked 
have been investigated earlier [2DH12]. Such a system 
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FIG. 2. The variation of Xc on the elastic modulus B of 
a phase separating binary mixture with elastic interactions 
for different values of migrant volume fraction (j>ao (0.05 (red 
dashed line), 0.10 (black solid line), 0.15 (blue dotted line) 
and 0.20 (green dash-dotted line)). Figure shows Xc increases 
with B (as ~ \/b) indicating that softer systems are more 
susceptible to phase separation and decreases with (pao for 
a fixed B. Inset shows phase diagram of polymer mixtures 
without elastic interactions. For X < Xc the system remains 
in one phase. Increasing x beyond Xc phase separation occurs. 
The coexistence (green dashed line) and spinodal (red solid 
line) curves demarcating phase boundaries is shown. 


shows microphase separation and is different from the 
functional proposed here (Eq. @)- 

Before discussing the surface segregation process we 
consider the bulk thermodynamic behaviour of the sys¬ 
tem described by Eq.(|^. This can be obtained easily by 
minimising the elastic EH free energy with respect to 4>. 
The minimisation procedure leads to a relation between 
(j) and Xi which for bulk concentration (poo corresponds 
to the binodal curve X = [1 ~ log(<^oo) — .^AA^]/[iVx(l — 
2(/)oo)]. It is interesting to note that x parameter does not 
depend on the elastic modulus B. The critical value of Xc 
above which the mixed phase is unstable, obtained from 
the relation d^Ffhe/dcp^ = 0, however increases with in¬ 
creasing B. This is shown in Fig.([^ with Xc ~ '•/b 
for different values of (poo- As shown in Fig. Xc de¬ 
creases with increasing (poo for a fixed B. This can be 
understood as follows. As the volume fraction of the 
migrant increases, the available free volume decreases 
and hence entropy decreases. Since a balance between 
entropic and enthalpic contributions dictates the equi¬ 
librium, a lower value of enthalpy (and hence lower x) 
is required to bring about the phase separation. With 
Xc increasing with B the single phase region of a rigid 
system is stable for larger values of x in comparison to 
polymer mixtures without elastic interactions. The phase 
behaviour of the binary polymer mixture without matrix 
elasticity is shown in the inset of Fig.([^. 


Surface segregation for elastic FH theory: The SB for¬ 
malism, outlined earlier can be used to compute the con¬ 
centration profile of the migrant (p{z) close to the in¬ 
terface for the phenomenological free energy functional 
described by Eq. (§. Fig.(§ shows migrant concentra¬ 
tion profiles for both systems, a symmetric binary poly¬ 
mer mixture having a bulk concentration (poo = 0.05 
with and without elastic interactions. The inset shows 
(p{z) as a function of depth z for different values of x 
for Nx = Nb = 10 in the absence of elasticity (ob¬ 
tained by minimising Eq.([^). For smaller values of x 
(—0.78 (red dashed line)) an approximate exponentially 
decaying profile is observed. As x increases, migrant 
concentration reaching the surface increases monoton- 
ically (x = 0.320,0.325) and beyond a critical value 
Xc = 0.327 a macroscopic wetting layer is observed. In 
contrast, when elastic interactions are included (main 
panel Fig.0), the migrant fraction for the same value 
of miscibility parameter x (0.320), obtained by integrat¬ 
ing the area under the curve (p{z) decreases monotoni- 
cally with increasing B. For lower values of the modu¬ 
lus, B = 0.1,0.108 a shoulder (reminiscent of a rounded 
wetting transition) is observed. For higher values of B 
(0.13, 0.3) an exponentially decaying profile is obtained 
suggesting elastic interactions severely inhibiting migra¬ 
tion. 

While physically intuitive and relatively straightfor¬ 
ward to implement the SB model has some disadvan¬ 
tages. First, the surface migrant fraction (pi is an addi¬ 
tional input and cannot be calculated from the model. In 
order to establish our main result, namely, elastic inter¬ 
actions inhibit surface migration as the matrix rigidity 
is increased, we employ a self-consistent field theoretic 
approach where this limitation does not exist. However, 
both the SB model and the SOFT framework suffers from 
the limitation that the bulk volume fraction (poo is held 
constant, no matter how much material flows to the sur¬ 
face. Modifications to the SB and SOFT framework that 
do not suffer from this drawback will be reported else¬ 
where |23) . 

Self-Consistent Field Theory: First introduced in con¬ 
text of polymers by Edwards [ 21 ], self-consistent field 
theory (SOFT) has been successfully employed to solve 
equilibrium behaviour of polymeric systems [25] . We em¬ 
ploy the SOFT formalism developed for end absorbed 
polymer brushes in polymer matrices HSIIIT] and adapt 
it to our situation. 

We consider two types of chains A, and B, with each 
component interacting via a mean-field potential on a 
lattice. The configuration of a polymer chain on the lat¬ 
tice can be visualised as a random walk with the variable 
representing a segment, akin to “time”. The distribu¬ 
tion function of a chain of a given species qk{z, t), (where 
k = A,B indicates the species) in the absence of interac¬ 
tion is Gaussian, satisfying a diffusion equation with the 
diffusion constant given by the squared radius of gyration 
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FIG. 3. Migrant concentration profiles 4>{z) for the SB model 
including the elasticity obtained by minimising Eq. for 
(x = 0.320 and Na = 10) and increasing B (0.1 (green solid 
line), 0.108 (red dash double dotted line), 0.13 (yellow dash 
dotted line), 0.3 (blue dashed line)). A wetting transition 
is not observed in this model. Inset, shows concentration 
profiles for SB model without elasticity for symmetric case 
Na = Nb = 10 for increasing x (—0.78 (red dashed line), 
0.320 (green dotted line), 0.325 (blue dash dotted line) and 
0.327 (black solid line with a break) indicating the formation 
of a macroscopic wetting layer. 

Rg = Na?/a. When the chains interact, qk{z,t) satisfies 
the modified diffusion equation 

= ^^V^gfc(z,t) - Wk{z)qk{z,t), (6) 

where a is the Kuhn length, N the degree of polymerisa¬ 
tion and Wk{z) is a mean field which takes into account 
the interactions between the monomer and its neighbour¬ 
ing chains of either species. Since we consider segrega¬ 
tion in one dimension the variable z in Eq. represents 
the distance from the surface, while “time” t indicates a 
particular segment of a chain. The mean field potential 
for the migrant species wa(z) is related to it’s chemical 
potential by 

wa{z) = - ksT log (j)A{z)) , (7) 

with the bare chemical potential ^j?a{z) = NAFfhe + 
Na4>b {dFfhe/d(j)A- dFfhe/d4>B)- For the matrix, we 
neglect the entropic contribution in Eq. 0. Thus the 
mean field potential wb{z) is the same as the bare chem¬ 
ical potential /r®, which has the same functional form as 
with A and B indices interchanged. The volume frac¬ 
tions (pk can be computed from the distribution function 
using 

g/3/ife(z) r-Nk 

4'k{z) = —-— / dt qk{z,t)qkiz,Nk-t), (8) 
-'''fc Jo 


where /ifc(z) = ^^{z)+Fs for the matrix and the migrant, 
with Fs (defined in units of l/ksT) being the surface 
free energy. Eq.([^-Eq.(|^ are solved self-consistently [Ml 
HZ] with the boundary conditions at the surface qA^z = 
0,t) = and qB{z = 0,t) = 0 for the migrant and 

the matrix respectively. The initial conditions qk{z,t = 
0) = 1 is used for both species. 

The concentration profile of the migrant as a function 
of distance from the surface (in units of Rg) obtained 
from the SOFT calculation is shown in Fig. 0 for a misci¬ 
bility parameter x = 0.22 and surface energy Fg = —2.0. 
The migrant polymer has a Kuhn length a = 1 and a 
degree of polymerisation Na = 10. As the elastic modu¬ 
lus of the matrix B is increased (from 0.001 to 0.11) the 
amount of material migrating to the surface decreases 
monotonically. In contrast to the SB model where pi is 
an additional input, ((^i = 1.0) in Fig.([^), it can be cal¬ 
culated within the SOFT framework. Fig. 0 shows pi 
decreasing monotonically with increasing B. The inset 
shows Fig. 0 the variation of the migrant concentration 
at the surface pi as function of B for different values of 
the surface energy Fg. The effect of elasticity on the mi¬ 
grant fraction pi is more pronounced for low values of B, 
(« 0 — 0.02). As expected pi decreases with increasing 
surface free energy Fg for a given value of B. For the 
elastic systems considered here, a wetting transition is 
not observed. A direct comparison between the param¬ 
eters in SB model and a variant of the SOFT method 
presented here [28] is currently underway. 



FIG. 4. Migrant concentration profiles p{z) for different elas¬ 
tic modulii B of the polymer matrix (0.001, (black solid line), 
0.01 (magenta dash double dotted line), 0.05 (green dotted 
line), 0.08 (yellow dash dotted line), and 0.11 (blue dashed 
line)). The amount of material flowing to the surface de¬ 
creases with increasing B. The dependence of the surface 
fraction pi as a function of B for different surface free energy 
Fg (—0.795 (black solid line), 0.0 (blue dotted line), 0.368 
(green dash dotted line), and 0.963 (magenta dashed line)) is 
shown in the inset. As expected the volume fraction decreases 
for system with higher Fg. 




























5 



FIG. 5. Cahn construction showing first order wetting transi¬ 
tion for the Flory Huggins free energy functional, Ffh (Eq. 0) 
(green dash dotted line). An intersection between Ffhi<l>i) 
and Fg{(j)i) (blue solid line) at three points intersects demar¬ 
cating areas Si and S 2 , such that Si > S 2 indicates a first 
order wetting transition. A similar graphical construction 
for the elastic Flory-Huggins functional Ffhe (Eq.([^)) with 
B — 0.17 (red dashed line) shows one intersection, indicating 
the absence of wetting transition. 

Cahn construction: A geometric way of demonstrating 
the absence of a wetting transition has been proposed by 
Cahn [HI IH] and applied in context of binary polymer 
mixtures |29j . A calculation of the surface migrant con¬ 
centration (j)i involves solving the equation 

F'M = VHhWiM, (9) 

where has the same meaning as Eq.(0, and F{(j)i) 
refers to the Ffh for SB model and Ffhe when elastic 
interactions are present. A graphical method of solving 
is shown in Fig. 0, plotting the surface F'^{(f>i) 
(blue solid line) for = —0.5, and g = 0.4, and bulk 
free energy contributions y/k{cj)i)F{cj)i) as a function of 
(f>i for a system with (red dashed line) and without (green 
dash dotted line) elastic interactions. In the absence of 
elasticity B = 0 the curves intersect at three points, de¬ 
marcating areas Si and S 2 such that 5'i > S 2 indicating 
a wetting transition. For a finite value of B (0.17 in 
Fig.(0) the wetting transition is absent PE]. 

Conclusion: In conclusion, we have analysed the role 
of matrix elasticity on the surface migration of small 
molecules in binary polymer mixtures using mean field 
and self-consistent field theories. We have shown that 
increasing the rigidity of the matrix leads to significant 
reduction of the migrant fraction on the surface. Fur¬ 
ther, by increasing the elastic modulus of the polymer 
matrix a wetting transition can be avoided for systems 
having miscibility parameters in the range that would 
otherwise have led to it. This provides a novel way of 
controlling surface migration in complex industrial for¬ 


mulations such as adhesives in hygiene products where 
surface migration leads to decreased product functional¬ 
ity. To the best of our knowledge, the only experimental 
system (despite significant differences) related to the the¬ 
ory presented here investigates segregation processes in 
polystyrene networks [20]. We hope that our theoretical 
work will prompt experimental studies in this direction. 

Acknowledgements: JK, TCBM, and BC thank Proc¬ 
ter and Gamble for financial support. SC and BC 
acknowledge funding via the ERC Marie Curie ITN 
MICSED grant under FP7 program. The authors thank 
Richard Thompson, Elise Sabattie, Mark Wilson and 
Jos Tasche from Durham University, Marc Hamm from 
Henkel, and Todd Mansfield, Gabriela Schaefer, Torsten 
Lindner and Mattias Schmidt from P&G for helpful dis¬ 
cussions. 


* buddhapriya.chakrabarti@durham.ac.uk 

[1] R. A. Jones and R. W. Richards, Polymers at Surfaces 
and Interfaces, (Cambridge University Press, 1999). 

[2] A. Budkowski, Adv. Poly. Sci. 148, 3 (1999). 

[3] Y. Lipatov, Prog. Polym. Sci. 27, 1721 (2002). 

[4] P. Lonchampt and R. W. Hartel, Eur. J. Lipid Sci. Tech- 
nol. 106, 241 (2004). 

[5] K. Bhunia, S. S. Sablani, J. Tang, and B. Rasco, Compr. 
Rev. Food Sci. F 12, 523 (2013). 

[6] M. Geoghegan and G. Kraush, Prog. Polym. Sci. 28, 261 
(2003). 

[7] K. Binder, Acta Polymer 46, 204 (1995). 

[8] J. Cahn, J. Chem. Phys. 66, 3667 (1977). 

[9] P.-G. de Gennes, Rev. Mod. Phys. 57, 827 (1985). 

[10] D.-K. Pan and W. M. Prest, J. Appl. Phys. 58, 2861 
(1985). 

[11] R. A. L. Jones et at, Europhysics Letters 12, 41 (1990). 

[12] I. Schmidt and K. Binder, J. Physique 46, 1631 (1985). 

[13] R. A. L. Jones et al, Phys. Rev. Lett. 62, 280 (1989). 

[14] U. Steiner et al., Science 258, 1126 (1992). 

[15] M. Rubinstein and R. Colby, Polymer Physics, 1st ed. 
(Oxford University Press, 2009). 

[16] P.-G. de Gennes, Scaling Concept in Polymer Physics, 
(Cornell University Press, Ithaca, New York, 1979). 

[17] R. A. L. Jones, Phys. Rev. E 47, 1437 (1993). 

[18] P. J. Flory and J. Rehner, J. Chem. Phys. 11, 521 (1943). 

[19] L. R. G. Treloar, The Physics of Rubber Elasticity, 3rd 
ed. (Clarendon Press, 2005). 

[20] P.-G. de Gennes, J. Physique 40, L69 (1979). 

[21] D. Read, M. Brereton, and T. C. B. McLeish, J.Phys II 
France 5, 1679 (1995). 

[22] N. Clarke, T. C. B. McLeish, S. Pavawongsak, and J. 
Higgins, Macromolecules 30, 4459 (1997). 

[23] J. Krawczyk, S. Croce, and B. Chakrabarti, (2015), un¬ 
published. 

[24] J. S. Edwards, Proc. Phys. Soc. 85, 613 (1965). 

[25] G. Fredrickson, Equilibrium Theory of Inhomogeneous 
Polymers, 1st ed. (Oxford University Press, 2006). 

[26] K. R. Shull and E. Kramer, Macromol. 23, 4769 (1990). 

[27] K. R. Shull, J. Chem. Phys. 94, 5723 (1991). 

[28] J. Genzer, A. Faldi, and R. J. Composto, Phys. Rev. E 















6 


50, 6 (1994). [30] M. Geoghegan et ai, J. Phys.: Condens. Matter 12, 5129 

[29] A. Budkowski, F. Scheffold, J. Klein, and L. J. Fetters, (2000). 

J. Chem. Phys. 106, 719 (1997). 



